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We study spin wave relaxation in quantum Hall ferromagnet regimes. Spin-orbit 
coupling is considered as a factor determining spin nonconservation, and external 
random potential as a cause of energy dissipation making spin-flip processes irre- 
versible. We compare this relaxation mechanism with other relaxation channels 
existing in a quantum Hall ferromagnet. 
PACS numbers 73.21.Fg, 73.43.Lp, 78.67.De 

1. Last years are characterized by growing interest in spin relaxation (SR) in low- 
dimension systems — first of all, in the relaxation in quantum dots studied within the 
projects aimed at development of a computer employing spin memory. Yet, the relaxation 
of an electron spin in lateral quantum dots manufactured on the basis of two-dimensional 
(2D) heterostructures, should be in many respects similar to the SR of electrons localized in 
the 2D layer in minima of a smooth random potential (SRP). In high magnetic fields this 
single-electron relaxation corresponds to the situation occurring at low Landau level (LL) 
filling: z/ <^ 1 or \v — 2n\ <^ 1 (n is an integer) .- 

The SR at different filing factors, has quite different nature representing in this 

case a many-electron process. In particular, in a quantum Hall ferromagnet (QHF), i.e. at 
u = 1,3,... or 1/ = 1/3,1/5,..., the SR reduces to the relaxation of lowest collective ex- 
citations, i.e. spin waves. ^''^ The SR observation would thereby be a good tool to study 
fundamental collective properties of a strongly correlated 2D electron gas (2DEG). However, 
in spite of much recent interest in the SR in a 2DEG, up to now only a handful of exper- 
iments relevant to the SR in a QHF were performed: these are indirect results based on 
the linewidth measurements in the electron spin resonance,- and a direct observation where 
the photoluminescence dynamics of spin-up and spin-down states was studied.- Meanwhile, 
availability of the new time-resolved technique of photon counting allows us to believe that 
new direct experiments on observation of excitations' relaxation in a 2DEG, in particular of 
the spin wave relaxation (SWR), will become available in the near future.- 

Theoretically the SWR in a QHF was studied in works 7,8. It is worth noting here that 
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the SWR represents actually not spin dephasing but the energy relaxation due to the spin-flip 
process. Indeed, any spin- flip means at least dissipation of the Zeeman energy ez = \g\fiBB 
[g —0.44 in a GaAs structure). The latter is a part of the spin- wave (spin exciton, SE) 
energy 

-Esw = ez + Sq-, (1) 

whe.e e, i= the SE collation energy depending on the 2D wave vecto. At wianee with 
the relaxation channel of Ref. |7| where electron-phonon interaction was considered as the 
mechanism making the relaxation irreversible, and contrary to the case of Ref. ^ where the 
irreversibility was provided by an inter-spin-exciton interaction mechanism, we now study 
smooth disorder field as the reason causing the energy transform. The SRP thereby deter- 
mines an alternative relaxation channel competing with the ones studied earlier. Another 
distinction of the present work from Refs. |7||8| consists in the study of not only the integer 
QHF (at z/ = 1, 3, ...) but also of the fractional one (z/ = 1/3, 1/5, ...) as well. At the same 
time we again consider the spin-orbit coupling (SO) as the cause mixing different spin states 
and therefore providing the spin nonconservation. Actually, various SWR channels coexist 
in parallel. We consider the total rate and find crossover regions of external parameters 
(magnetic field, temperature, etc.) where one relaxation channel ceases to be dominant and 
changes into another. 

The SR channel due to SRP was already considered in the integer quantum Hall ferro- 
magnetic case.^'^ However, studied in these works instead of the SWR was a specific SR 
when initially the total macroscopic spin S of the system as a whole is turned away from 
the equilibrium direction parallel to B. (Relaxation of this Goldstone mode microscopically 
reduces to annihilation processes of the so-called zero SEs, having exactly zero momenta.) 
Contrary to this case, the spin perturbation determined by excitation of the spin waves 
(non-zero SEs) represents an initial deviation where AS = ASz-, so that S is kept parallel to 
B and the total symmetry of system remains unchanged. 

Concerning the origin of SRP, one should note that it has in the 2D layer the "direct" 
component and the effective one. The former is the SRP determined by charged donors 
located outside the spacer. The latter is essential in some kinds of quantum wells, being 
determined by spatial fluctuations (in the plane of the layer) of quantum well width. These 
fluctuations lead to fluctuations of the size-quantization energy and may be presented as an 
SRP term in the single electron Hamiltonian. Both SRP components have approximately 
the same amplitude A~10K and correlation length A~30 — 50 nm. 

2. The total Hamiltonian has form iftot = Xli-^i"'^ + -^int; where j enumerates electrons. 
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i^int is the e-e interaction, and the single-electron operator is 

= h\'/2m: - ezaj2 + Hso + ^(r) . (2) 

In this equation ip{r) is the SRP field; the SO Hamiltonian is specified for the (001) GaAs 
plane, 

Hso = a (q X (t)^+/3 {qyCTy-q^ax) , (3) 

presenting a combination of the Rashba term and the crystalline anisotropy term — (q = 
— zV + eA/ch is a 2D operator, cr^.,^,^ are the Pauli matrices). If the SRP is assumed to be 
Gaussian, then it is defined by the correlator K{r) = ((y9(r)(/?(0)). By choosing {^p{r)) = 0, 
in terms of the correlation length A and the LL width A the correlator is 

K{r) = exp (-rVA') • (4) 

We first find the bare single-electron basis diagonalizing the Hamiltonian (2) without the 
SRP field. To within the leading order in the Hso terms we obtain 

X^J = ( 'I'^v \ 

/ ^ ^ \ 

\^ ^np j 

Here ipnp is the electron wave function in the Landau gauge, n is the number of the half-filled 
LL in the odd-integer quantum Hall regime, i.e. in the v = 2n+l case. Otherwise, if 1, we 
set n = 0. u and v are small dimensionless parameters: u = (3\/2/lBhoJc and v = a\/2/lBhuJc 
{(jjc and Ib are the cyclotron frequency and the magnetic length, respectively). The single- 
electron states thus cease to be purely spin states but acquire a chirality a or h. The spin 
fiip corresponds thereby to the a ^ b process now. 



By analogy with previous works-'^*^'^ (see also Ref. I 111 ) we define the SE creation operator 



where ap and bp are the Fermi annihilation operators corresponding to states (5), is 
the LL degeneracy number. In Eq. (6) and everywhere below we measure wave vector q 
in the 1/Ib units. If the ratio = {ae^ / hIb) /^ujc is considered to be small (a < 1 is the 
averaged formfactor which appears due to finiteness of the layer thickness), and the SRP 
and SO terms in Eq. (2) are ignored, then the operator (6) acting on the ground state in 
the odd-integer quantum Hall regime yields the eigen state of the total Hamiltonian: namely, 

[Htou QLq]|0) = (ez+^^g)Qlfeq|0), wherc |0) = I T, T, ••• T )• This basic property of the exciton 
state, QlbqiO)' is asymptotically exact one to the first order in r^. 
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Now consider corrections arising due to the Hso terms. When presented in terms of basis 
states (5), spin operators J \l/"''S^\l/(i^r and J '^^Sz'^d^r [where \l/ = ^p(ap\l/pa-|-6p\l/pf,)] preserve 
invariant form up to the second order in u and v. However, the interaction Hamiltonian 
-f^int = ^ /(iri(ir2 \[^^(r2)\I/^(ri)f/(ri — r2)\[^(ri)\&(r2) acquires proportional to u and v terms 
which correspond to creation and annihilation of SEs in the system. It is exactly these terms 
that lead to the "coalescence" channel of the SWR.- In the present work we study another 
relaxation channel. Therefore, neglecting this SO corrections to i^int, we focus on the SRP 
term. Calculating J \E'"'"(y9(r)\E'(i^r, we get the terms responsible for a spin- flip: 

^ = Nl^%J2^{<l) {^uq+ - vq^) + H.c. (7) 
q 

(it is assumed here that g^l). '^(q) is the Fourier component [i.e. = X]q^(^)^*''1; ^'^d 

?± = =F^(gx±^gj/)/y2- 

At variance with integer QHF, the use of the excitonic basis Q\ij^^) presents only a 
model approach in the case of fractional quantum Hall regime. Generally, spin-flip excitations 
within the same Landau level might be many-particle rather than two-particle excitations 
at fractional filling because the same change of the spin numbers 5S = 6Sz = — 1 may be 
achieved with participation of arbitrary number of intra-spin-sublevel excitations (charge- 
density waves). These waves are generated by the operator Al^ = N^^^'^ Ql^^ acting on the 

ground state |0) = | t, •• T,-- T)-^^ It is trivial in the case of integer z/ (^q|0) = 5q,o|0)); 
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however, states of the Qlbq^-^qa-^qa"-!'-') ^^P^ might constitute a basis set if one studies 
spin-flip at fractional u. On the other hand, a comprehensive phenomenological analysis-ii^ 
suggests that even the spin-flip basis reduced to single-mode (single-exciton) states would be 
quite appropriate, at least for lowest-energy excitations in the case of fractional QHF. This 
single-mode approach is indirectly substantiated by the fact that the charge-density wave has 
a Coulomb gap'^ which is well larger than the Zeeman gap ez- Hence for a fractional QHF, 
just as in Ref. 3!, we will consider the only state Q|j5q|0) to describe the spin-flip excitation. 
The commutation algebra for operators Q^^^q, A]^, and ^B^,, = N^^^'^Qlf^^,, is certainly the 
same as for integer filling,-'^'^. However, a difference arises in the calculation of expectation 
(0|^q^q,|0) which is needful for the following. This value is simply 5q,o5q',o at integer filling, 
but at z/< 1 it is expressed in terms of the two-particle correlation function g{r) calculated 
for the ground state: 

{0\A^A[,\0) = ^ [27rugiq)e'^'/' + l\ <5q, q . (8) 

Here 'g{q) = f g{r)e~'^'^^d'^r is the Fourier component. Function g{r) is well known, e.g., in 
the case of Laughlin's state.— If the ground state is presented in terms of the Hartree-Fock 
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model, we get the expression 2'7ig= (^N^5g^o — e~'^'^^^ which does not depend on z/. Besides, 
at odd-integer filhng factors this Hartree-Fock expression becomes Fourier component of the 
exact correlation function. In the latter case one should also make the substitution u v—2n 
in Eq. (8), i.e. formally set u = l there. 

3. The operator (7) obviously does not conserve the number of SEs. However, if the 
SWR is governed by this operator, the corresponding problem can not be solved in terms 
of a single-exciton study. Indeed, the SE interaction with the SRP incorporates the energy 
f4-sRP ~ Q'^B^/A (the SE possesses the dipole momentum e/B[qxi])^. The SE momentum 
is estimated from the condition Sq<T, and we therefore find that U^sbp ^ ^z, T. Due 
to this inequality, the energy of annihilating exciton can not be transformed to anywhere. 
By analogy with Ref. sl, we study a coalescence process where initial double-exciton state 
\i) = SlbqiSlfeqalO) transforms to final single-exciton state |/) = Qlbq/|0) having the combined 
energy: 

ez + Sg> = 2ez + Sg, + (9) 

(c.f. also the Auger magnetoplasma relaxation considered in Ref. \\\ . At the same time, 
contrary to Ref. la, there is no momentum conservation in this SWR channel. Thus the phase 
volume where the Xq^+Xqj —^Xqi transition is possible turns out to be much larger than 
that in the coalescence process of Ref. ^. This transition is governed by the Fermi golden 
rule probability: Wfi = {2TT/h)\/iAfi\'^6{Ef — Ei), and our immediate task is to calculate the 
matrix element Ai fi = u^'^^'^{f\(p\i). (The factor z/^'^/^ appears due to the normalization since 
norms of the \i) and |/) states are u"^ and z/, respectively.) 

We perform the calculation for relevant values of momenta qi,q2,q' <^ 1 which satisfy 
the conditions Sq^,Sq^<T < IK. (These inequalities correspond to qi,q2,q' '^^/Ib in usual 
dimensional units). By employing exciton-operators' commutation rules - and evident iden- 
tities Qabq|0) = i3q|0)=0 and (0|^q|0) = z/, we obtain with the help of Eqs. (7)-(8) that 



-^/i(qi,q2,q) = — 



2 



^(^(q)(mg+-wg_)5q,4^2,qi^/. (10) 



Besides, within our approximation, g(q)e'^^^'^ should be replaced with g(q)e'^^^'^ . The 
latter quantity is equal to —1/2tt in the Hartree-Fock approach or — l/27rz/ when calculated 
in the case of Laughlin's ground state describing the fractional QHF. So, for = 1, 1/3, 1/5, 
replacing the terms in square brackets with — l/vrz/, we obtain a simple result: 



l-M/.(qi,q2,q')|' = 47rir(g) '^'^'''^'''^ 



(11) 

CFCii+q2-q' 
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It is used that the squared modulus of ^(q) may be expressed in terms of Fourier component 
of the correlator (4): \Tp{(\)\^ = 2TiK{q) /N^f,. In the Hartree-Fock model the expression (11) 
should be multiplied by v"^] therefore the calculated relaxation rate would be by a factor of 
v"^ slower. Notice also that if z/ = 3,5, one should formally set z/=l in Eqs. (10) and (11). 

The SWR rate is defined as the difference between the fluxes of annihilating and cre- 
ated SEs. We assume that the thermodynamic equilibrium in the system of spin waves is 
established much faster than the spin-flip processes occur so that the rate is 

^= ^ 5Z X '-^^^^^i' ^(^1 + ^2-^0 [nm^ (1 + n') - n' (1 + m) (1 + ^2)] . (12) 

qi. q2, q' 

The notations used here are Ei = ez+Sq^, ni = n{Ei) {i = 1, 2) and E' = ez+Sq', n' = n{E'), 
where the Bose distribution function is n{E) = l/(e^^^^^'^ — 1). The rate R is completely 
determined by Eqs. (11)-(12) and is a function of parameters B, T, and of the total number 
of SWs in the system: A^x= Sq^^l^z+^g)- In our case, when temperature is rather low, we 
can certainly use quadratic approximation for the "kinetic" exciton energy: EqK, /2M^. 
Chemical potential yU is determined by the ratio of the exciton number and the total spin: 
N^ifi) = — S . Calculating the quantity Nx'^=N^ , one obtains the equilibrium 

number of excitons. We will find the rate at the final stage of the relaxation process where 
Nyr—N^^^ <^Ni^\ So, by employing the quadratic approximation for the SE kinetic energy, and 
changing in Eqs. (11)-(12) from summations to integrations we obtain R= (^N^ — Ni^'^^ /^sTp, 
where 

l/rs.p= ^^'+;J^^- (e--/^-e--/^) Fskp(A^M.T/;|, ez/T) . (13) 

Here FsRp(a, P) is a dimensionless function arising as a result of integrations over qi and q2 
and averaging over angles = qi A q' and ^2 = q2 A q': 

e~^~^ dxdy 



'SRP 



>, (3) 




oio (1-e— /5)(l-e~2^/3)(i_g-..-^2/3) 
X dOi I de2r{x,y,ei,e2)exp[-ar{x,y,ei,e2)], 



where r(x, y, 61,62) =x+y+(3/2 — y/x+y+f3 {^/x cos 61 + ^/y cos 62) + y/xy cos {61 — 62). 



4. Now we calculate the numerical value of l/rgrp at typical SRP parameters and com- 
pare it with inverse relaxation times l/r^^e ai^d 1/Tph governed by the inter-SEs' interaction 
mechanism- and the SE-acoustic-phonon coupling.- We carry out this analysis for the 1/ = ! 
QHF assuming that A = 10 K and A = 40 nm. The Zeeman splitting a.t g = —0.44 is 
ez = 0.2955 K {B is everywhere in Teslas), and the combination of SO parameters is esti- 
mated as u^ + v"^ = 10^^ /B. The SE mass might be calculated theoretically by using 
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general expressions for Eq.-^ Yet, the result depends on specific formfactor inherent in a 
given heterostructure due to finite thickness and it is therefore more convenient to extract 
Mx immediately from experiments. According to recent data available for currently used 
wide quantum wells,— we estimate that 1/Mx = 9.24-\/5 K. Using Eq. (13), we thus 
calculate l/xsi-p as a function of temperature T at given field B. The results are presented 
in Fig. 1 by dash curves. The dot and dash-dot curves correspond to the 1/Te-e and l/rph 
values given by formulas— 

= l^u' + v')T (e--/^-e"^-/^) F,_,(ez/T) , (14) 



where 



and 



where 



F^,(/?) = jj 



dxdy {x + y + (3)e 



{xy-P^/Af'^ (1 - e-^) (1 - e-^) (1 - e-^M) 

xj/>/32/4 



-1 



(15) 







(See Ref. |7|; the used material parameters characterizing the electron-phonon coupling are 
c, = 5.14-105cm/s, = 0.8- 10^12 g-i^ rp = 35- lO^^^ g-i^ and po = 2.52-106cm-^ both kinds 
of e-'ph interaction, deformation and polarization ones, are taken into account.) 

It is seen from Fig. 1 that the SRP relaxation channel actually competes with other 
mechanisms in the experimentally relevant range of parameters: namely, at fields B < 5 
and temperatures T ~ 0.3 — 0.5 K. We have indicated above that the basic advantage of the 
SRP channel, as compared to the e-e one, consist in the absence of momentum conservation 
in the coalescence process. On the other hand, the SRP mechanisms is also determined by 
effective SE-SE collisions. Therefore the inverse relaxation time is proportional to the SE 
concentration and drops exponentially as ~exp (— ez/T) with vanishing T [rather than as ~ 
exp (— 2ez/T) which occurs for the e-e mechanism due to the SEs momentum conservation!]. 
The phonon mechanism of SWR dominates at low temperatures due to its weak temperature 
dependence (~T), in spite of small value of the electron-phonon coupling constant in GaAs. 
The dependence on the filling factor in the case of integer QHF is only determined by the 
SE mass because u in Eq. (13) is formally set equal to unit. For fractional QHF there 
are both direct and indirect (through the mass Mx) dependences on u. 

Finally we calculate the combined inverse relaxation time determined by the SO interac- 
tion: 

l/^ot = l/rsrp + l/T-e-e + l/7-pft (16) 
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The result is presented by solid curves in Fig. 1. It is worth mentioning that it demonstrates a 
good agreement with the measured value Ttot — 10 ns ofRef. \d when calculated for parameters 
B and T corresponding to the experiment. 
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1'' The corresponding formula for inverse relaxation time in Ref. 
present the corrected result in Eq. (14). 
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FIG. 1: Inverse SWR times against T calculated by using formulas (13)-(15) at B = 3, 5, 10 T. 
Specific material parameters are given in the text. Dash, dot, and dash-dot lines are for l/rgrp, 
and 1/Tph, respectively. Solid lines present the result of calculation of the combined inverse 
time (16). 



